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Abstract. Although hydrogen cyanide has become quite a common molecular 
tracing species for a variety of astrophysical sources, it, however, exhibits dra- 
matic non-LTE behaviour in its hyperfine line structure. Individual hyperfine 
components can be strongly boosted or suppressed. If these so-called hyper- 
fine line anomalies are present in the HCN rotational spectra towards low or 
high mass cores, this will affect the interpretation of various physical proper- 
ties such as the line opacity, ta, and excitation temperature, Tex, in the case 
of low mass objects and infall velocities, Vin, in the case of their higher mass 
counterparts. This is as a consequence of the direct affects that anomalies have 
on the underlying line shape, be it with the line structural width or through 
the inferred line strength. This work involves the first observational investiga- 
tion of these anomalies in two HCN rotational transitions, J=l— >0 and J=3^'2 
towards both low mass starless cores and high mass protostellar objects. The 
degree of anomaly in these two rotational transitions is considered by com- 
puting the ratios of neighboring hyperfine lines in individual spectra. Results 
indicate some degree of anomaly is present in all cores considered in our survey, 
the most likely cause being line overlap effects among hyperfine components in 
higher rotational transitions. 

Key words: Starless cores - Hyperfine anomalies - ISM: kinematics and dy- 
namics - Molecular Clouds - parameter space - dense tracers. 

1. INTRODUCTION 

To be a good probing species for dynamical conditions that lead to the onset of 
star formation, a molecule must be abundant enough so that it can be detected as 
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well as excitable in higher density regions. The purpose of the latter is such that 
it can be used to trace the deep dense inner regions rather than the sufficiently 
lower density outer regions. Being a nitrogen-bearing species, HCN fulfills these 
properties and having first been discovered in the cosmos by Snyder & Buhl (1979), 
it has since been used routinely as an infall tracer in low mass starless cores (Tafalla 
et al. 2006; Sohn et al. 2007) and also in high mass star formation regions (Wu 
et al. 2005). 

HCN has a hyperfinc structure due to the nuclear quadrupole moment of the 
^^N. This extends the use of a particular species since the optical depth and 
self-absorption could be measured by examining individual hyperfine lines while 
a measure of the excitation temperature may be gleaned from an examination of 
different rotational transitions. 

Low mass J=l-0 
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Fig. 1. Top (left): Optically thin HCN J=l— >0 for a low mass quiescent cloud core. 
Top (right): TRAO observation of low mass protostellar core, L1197. Bottom [left): 
Optically thin HCN J=l— 7>0 for high mass turbulent cloud core. Bottom (right): Mopra 
observation of high mass core in the G333-6.02 molecular cloud. 



Firstly identified by Gottlieb et al. (1975), the hyperfine anomalies present 
in HCN J=1^0 were confirmed towards the low mass starless core TMC-1 by 
Walmslcy et al. (1982) who demonstrated that the hyperfine components in the 
core spectrum were present in ratios that were not in thermal equilibrium with one 
another. Such anomalies are seen to be present towards mostly nitrogen-bearing 
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Fig. 2. Left: Optically thin HCN J=3— >-2 compared with Right: a JCMT observation 
of L1622A-2. 

species, i.e. molecules possessing the ^"^N-isotope, where the individual hyper- 
fine lines can be seen to be boosted or suppressed far beyond LTE or saturation 
conditions. 

2. HCN AS AN OBSERVING TOOL 

The lower rotational transitions of HCN arc excellent tracers of dense molecular 
gas in star-forming clouds. The relation Ucrit oc iJ,^i'j^j_i (Papadopolous 2007) 
for critical densities of rotational transitions in which i/j^j-i is the frequency of 
an optically thin line, allow HCN transition lines to trace approximately 100-500 
times denser gas than corresponding CO transitions. 

HCN is less prone to freeze-out compared to other abundant H2 mass tracers 
such as HCO+ or C^''0 (Freed & Mangum 2005), and therefore remains abundant 
in the gas phase in the cold central regions of star-forming cores. The hyperfine 
structure of HCN is due solely to the non- vanishing electric quadrupole moment 
of the end ^^N nucleus, leading to line splitting of several MHz, resulting in an 
easily identified split line structure. Hyperfine line ratios prove an indispensable 
tool in the constraining of optical depth estimations of a given region in the cloud. 

As mentioned in the foregoing section, anomalies scupper the usefulness of 
HCN in tracing the dynamics of interstellar regions. This is due principally to 
the variation in relative intensity among hyperfine components. Under LTE or 
optically thin conditions, the relative weightings of the hyperfine components in 
the lower rotational transition are of the form 1:5:3. For HCN J=3— >2, four of 
the six hyperfine lines are not spectrally resolved and the spectrum takes the 
appearance of one strong central component and two satellite lines, with a ratio 
of 1:25:1 between them. 

In Figure 1, the line marked "A" (JF=10— >-01) in the low-mass spectrum is 
boosted far beyond its expected strength. As for the high-mass spectrum in the 
bottom panel, given that the hyperfine components are broader on account of 
the increased turbulence, the component marked "A" is again seen to be boosted 
as well as broadened relative to the other components. This strongly suggests 
that line overlap effects at higher energy levels are significant. For our analysis 
of the HCN J=3— >2 hyperfine spectrum, wc defined the leftward component as 
the AF=0~ hyperfine branch, the rightward component as the AF=0"'" hyperfine 
branch and the central four components that remain unresolved towards even the 
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most quiescent low mass cores as the AF=1 hyperfine branch. This nomencla- 
ture is in keeping with the relative changes in hyperfine quantum number for the 
components involved (except the JF=32^>23, where AF=-1, which is part of the 
central group of 4 components) as well as the relative velocity shifts of the two 
AF=0 components, (-)-leftward versus (+)-rightward. In relation to Figure 2, the 
1:25:1 proportion of intensity amongst the hyperfine branches in optically thin con- 
ditions is seen to be unmatched towards the low-mass source with the component 
marked "B" (JF=33^23) appearing strongly anomalous. 

3. OBSERVATIONAL DATA 

In order to quantify the presence of HCN hyperfine anomalies in low mass star- 
less cores, observations of both the HCN J=l— >0 and J=3— >-2 lines for a selection 
of low-mass starless sources, mainly towards the Tuarus- Auriga and Ophiuchus 
molecular clouds, were collected. The sources were primarily chosen from a com- 
prehensive HCN J=l— >0 survey in HCN and HNC, catalogued by Sohn et al. 
(2007) using the Tacduk Radio Astronomical Observatory (TRAO). The low- mass 
HCN data forming part of this work is comprised of both the TRAO data ob- 
served by Sohn ct al. (2007) of 65 sources observed in J=l— >-0 and 29 of these 
same sources observed in J=3— >-2. 

The cores observed in J=3^2 were selected on the basis of how bright the 
central F=2— >1 hyperfine component is in the HCN J=l— >0 rotational transition 
since it was anticipated that if this particular hyperfine transition was strong then 
the probability for a detection in the J=3^'2 rotational transition of HCN was im- 
proved. Cores that were also clearly anomalous in the HCN J=l— >^0 line were fur- 
ther prioritized. These observations were carried out at the James Clerk Maxwell 
Telescope (JCMT) between September 2007 and July 2008 in band 5 (T225 > 0.2) 
conditions. Single-point observations were obtained for 30 sources in position- 
switching mode with a predefined off-position of 500" in an arbitrary direction. 
The ACSIS digital autocorrelation spectrometer with a bandwidth of 250-MHz 
was used providing a velocity resolution of ~ 0.034kms~^. The receiver noise 
temperature (DSB mode) was 510-850K. The telescope main-beam efficiency was 
0.69 and the half-power beamwidth (HPBW) was ~20" in the range 225-285GHz. 
The reduced J=3— >2 data of each low-mass source along with its corresponding 
J=l— >-0 spectrum will be published in Loughnane et al. (2011). 

For the high mass HCN J=l— fO data, individual sources in the G333 molecular 
cloud associated with the RCW 106 HII region were observed (Lo et al. 2009). 
This giant molecular cloud complex spans an ~ 0.7 deg^ region and is ~ 3.6 kpc 
away. The data analysed as part of this work was collected with the 8-GHz wide 
band Mopra spectrometer (MOPS) centred at 87-GHz. The FWHM beamsize is 
^ 36" (Ladd et al. 2005). Full details of the observational setup is available in 
Lo et al. (2009). We used HNC J=l— >-0 data to locate single sources of emission 
in this data since this species is co-spatial to HCN but is less abundant. Maximal 
positions of intensity in the HCN J=1^0 line that correspond to a single peak 
in the HNC J=l— !>0 line (actually exhibits an unresolved hyperfine structure) 
represent our seven young massive protostellar cores as indicated in the integrated 
intensity zeroth moment maps in Figure 3. The reduced core spectra for these 
seven sources are displayed in Figure 5 of Loughnane et al. (2011). 
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Fig. 3. C^**0 J=1^>0 (left) and HCO+ J=l->0 (right) integrated intensity 
(zeroth-moment) maps of the G333 giant molecular cloud. The crosses mark the 
positions of the seven anomalous HCN emission sources presented in this work. 



4. RESULTS & CONCLUSIONS 

From the high mass HCN J=l— >-0 spectrum in Figure 2 (Mopra data), it is 
clear that the three hyperfine components are merged to a degree. In order to 
analyse the degree of anomaly in these high mass sources, it is therefore necessary 
to come up with a method of dealing with the integrated emission that falls in 
each overlapped region. Therefore for each high mass spectrum, the individual 
components were approximated as gaussian profiles with associated function given 
by Equation [T] below 



y{x) = Aexp —0.5 



(1) 



Rearranging Equation [TJ we were able to compute the velocity Vmid, defined 
as the velocity position between two hyperfine components where the height (or 
strength) of the integrated emission in the common region is at its greatest. The 
velocity, Vmid, is the solution to the quadratic equation given by Equation[5]below. 
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Using this velocity and a calculation of the integrated emission in each of the 
estimated gaussian profiles with similar linewidths and line strengths as the actual 
observed hyperfine components, the common intensity regions between two high 
mass hyperfine lines could be calculated. The details of this calculation are left 
to the appendices of Loughnane et al. (2011). The common intensity in each 
case was apportioned to each of the overlapping hyperfine lines by way of their 
non-overlapping intensities. In this way, an integrated emission calculation was 
computed for each hyperfine line belonging to an individual high mass core. 
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Fig. 4. HCN J=l— >-0 hyperfine ratios, R12 plotted against R02, for (a) 64 low-mass 
starless sources located in the Taurus-Auriga and Ophiuchus molecular clouds and (b) 
15 high-mass objects considered in this work, made up from seven high-mass protostcUar 
cores located in G333-6.02 and the eight sources in figure 1 of Pirogov (1999). In (c), 
both sets of object's ratios are plotted in tandem. In each of the three panels, the dot- 
dashed curve marked LTE represents the line along which cores displaying HCN hyperfine 
spectral components obeying the Boltzmann equations of thermodynamic equilibrium 
would lie, in that part of the cloud; hence local-thermodynamic equilibrium (LTE). 



For each of the seven sources observed in the G333 molecular cloud, we com- 
puted the relative intensity ratios R02; the ratio of the intensity of the F=0— )■! 
hyperfine component to the ratio of the intensity of the F=2— >1 hyperfine com- 
ponent, and R12; the ratio of the intensity of the F=l^'l component to the ratio 
of the intensity of the F=2— >1 hyperfine component. These sets of intensities for 
all seven sources were plotted against one another along with the same ratios for 
eight sources considered as part of Pirogov (1999). The resulting plot is displayed 
in the top right-hand panel of Figmc 4 (Fig. 4(b)) below. The same analysis 
was carried out for 65 of the low-mass sources considered in Sohn et al. (2007), 
with the resulting plot displayed in the top left-hand panel of Figure 4 (Fig. 4(a)) 
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below. Both sets of data points are included in tandem in the bottom panel (Fig. 

4(c)). 

What this analysis demonstrates is that high mass cores arc more likely to ex- 
hibit intensity ratios close to the 1:5:3 formation common to optically thin regions. 
In Figure 4(b), the high mass cores arc more tightly packed than the low mass 
objects in Figure 4(a), which show degrees of anomaly that are quite varied. In 
each of these panels, the; dashed curve represents the progression of LTE ratios for 
different optical thicknesses. 
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Fig. 5. Hyperfine ratios R12 {top) and R02 {bottom) versus N2H+ and N2D"'" column 
densities (both J=1^0 transition). The correlation coefficient r is recorded in the top 
left corner. 

Taking the sets of ratios computed above, individually, we set about comparing 
the degree of anomaly for low-mass starless cores with a physical parameter of the 
core. The columnar densities of both the N2H+ and N2D+ J=l— >0 transitions 
were taken from Crapsi et al. (2005) and plotted against the hyperfine ratios for 
the corresponding cores in our observational set. The results of these comparisons 
are shown in the plots of Figure 5 where the linear correlation coefficient for each 
plot is also given. From this analysis, although a linear correlation may be far too 
simplified, it is clear that there is a minor correlation between degree of anomaly 
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and increasing density of the central source. How this correlation impacts the 
formation of the hyperfine anomalies is unclear and since in mappings of several 
of our cores in HCN J=3— >2, there appears to be a sustaining of the appearance 
of the anomalous profile where the density would be lower, it is clear that there 
are other physical factors contributing to the anomalous lineshape. 
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